Photoionization of a strongly polarizable target 



A V Korol^'^ and A V Solov'yov^ 

^ Frankfurt Institute for Advanced Studies, Johann Wolfgang Goethe-Universitat, 
Ruth-Moufang-Str. 1, 60438 Frankfurt am Main, Germany 
^ Department of Physics, St. Petersburg State Maritime Technical University, 
Leninskii prospect 101, St. Petersburg 198262, Russia 

E-mail: korol@fias.uni-frankfurt.de, solovyovQfias.uni-frankfurt.de 

Abstract. We demonstrate that the angular distribution of photoelectrons from 
a strongly polarizable target exposed to a laser field can deviate noticeably from the 
prediction of conventional theory. Even within the dipole-photon approximation the 
profile of distribution is modified due to the action of the field of alternating dipole 
moment induced at the residue by the laser field. This effect, being quite sensitive to the 
dynamic polarizability of the residue and to its geometry, depends also on the intensity 
and frequency of the laser field. Numerical results, presented for sodium cluster anions, 
demonstrate that dramatic changes to the profile occur for the photon energies in 
vicinities of the plasmon resonances, where the effect is enhanced due to the increase 
in the residue polarizability. Strong modifications of the characteristics of a single- 
photon ionization process can be achieved by applying laser fields of comparatively 
low intensities Iq ^ 10^° - 10" W/cm^. 



1. Introduction 

In this paper we demonstrate that the profile of angular distribution of electrons emitted 
in the process of a single-photon ionization from a strongly polarizable target (in 
particular, a metal cluster anion) exposed to a laser field can be noticeably modified 
due to the action of the field U{r,t) of alternating dipole moment D(t) induced at the 
residue by the laser field. The field U (r, t) acts on the escaping electron and, being 
dependent on the polarization vector of the laser field, brings additional dependence of 
the cross section on the emission angle. The degree of this dependence is determined 
by the magnitude of D(t) which, in turn, is proportional to the residue's dynamic 
polarizability. Therefore, for a target with large polarizability the effect can be very 
pronounced even for laser fields of low intensities. 

The effect of the laser-field-induced dipole moment in atomic photoionization 
process was discussed in [1] (where one finds also the references to the earlier works 
by the authors) and, more recently, in [2]. Both papers were devoted to the study of 
multiphoton detachment of electrons from atomic negative ions. In [1] the process was 
analyzed by means of the perturbation theory. Numerical analysis was carried out for 
halogen and alkali metal anions, and it demonstrated, that the induced dipole moment 
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influences the magnitude of total cross section of the many-photon detachment. The 
paper [2] dealt with the process of a strong-field many-photon detachment of an atomic 
anion. The non-perturbative adiabatic theory, developed in [3] , was modified to account 
for the field of the induced dipole moment. The numerical calculations, performed for 
Rb~ and Cs~, showed the increase of the detachment rate and essential modification of 
the angular distribution. 

In our paper we investigate the influence of the induced dipole moment on the 
angular and spectral distributions of photoelectrons in the process of single-photon 
ionization (or detachment) occurring in the weak-field regime. Numerical results are 
presented for sodium cluster anions Na^ with = 10^ . . . 10^. The corresponding 
residue clusters are characterized by large values of the dynamic dipole polarizability 
a{Lj), which in the vicinity of plasmon resonances, i.e. at a; = 2 ... 3 eV, can be as large 
as 10^ . . . 10^ a.u., which by far exceed the maximum values of a{uj) of individual atoms 
(e.g., a{Lj) ^ 400 a.u. for Cs at the photon energy 0.47 eV [2]). 

We demonstrate that dramatic changes to the profile of angular distribution occur in 
vicinities of the plasmon resonances and for comparatively low intensities ~ 10^*^ — 10^^ 
W/cm^ of the laser field. These are exactly the ranges used in recent experimental 
studies [4, 5] of the photo detachment from various metal cluster anions carried out by 
means of angle resolved photoelectron spectroscopy. Therefore, the experimental test of 
the predicted effects seems to be feasible. 

The atomic system of units, e = m = h= lis used throughout the paper. 



2. Theoretical framework 
2.1. Setting the problem 

Consider the process of photoionization of a strongly polarizable target (a cluster, a 
fuUerene, a complex molecule, etc) by spatially homogenous linearly polarized laser 
field 

E(t) = Eo cos ut = ^ (Eoe"'"* + Eoe'"*) . (1) 

We assume, that the field wavelength A greatly exceeds the characteristic size R of the 
target, R. Thus, the process can be treated within the dipole-photon limit. 

The laser field ionizes the target and polarizes the residue. Therefore, in addition to 
other fields (the static potential of the residue, the long-range polarizational potential) 
the escaping photoelectron feels the field due to the oscillating dipole moment D(i) 
induced at the residue: 

(2) 

To simplify the consideration we assume, that the residue is spherically symmetric, 
so that its tensor of polarizability contains only a scalar part which is notated as a (a;). 
Then, the moment D(t) reads 

D(t) = de-'"* + d*e'"*, (3) 
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where d = d{uj) = a{uj) Eo/2. The field of the induced dipole can be written as follows: 

W(r, t) = U{r) e^''^* + U*{r) e'-^*, (4) 

with U{t) = dr/r^. This field is dependent on the direction of Eq, therefore, the angular 
distribution of photoelectrons acquires additional dependence on the polarization of the 
laser field. The degree of this dependence is determined by the magnitude of D(t), which 
is proportional to the residue dynamic dipole polarizability. Hence, for a target with 
large polarizability the modification of the angular distribution due to the induced dipole 
moment can become noticeable even for a laser field of a comparatively low intensity. 

2.2. The perturbation series 

Let us clearly state the conditions which we impose on the parameters of the laser field. 
These include the following inequalities: 

(I)Eo<Eint, (II) £>^, (III)c^>|£i|. (5) 

The first condition implies that the laser field strength is much smaller than the strength 
-Eint of the internal field in the target. This condition allows one to ignore the laser field 
when constructing the wavefunction of the bound electron with the binding energy 
< 0. Numerical results presented in section 3 refer to the laser field strength 
Eq ~ 10-3 a.u. ~ 10^ V/cm. 

The second condition states that the energy e of the escaping electron is large 
compared to the ponderomotive shift due to the wiggling in the laser field. Owing to 
this condition we do not "dress" the photoelectron wavefunction with the laser field. 

The last inequality in (5) indicates that the photon energy exceeds the ionization 
threshold of the bound electron. Therefore, the process of a single-photon ionization is 
allowed for which the energy conservation law reads 

e = Ei + uj . (6) 

Combining the conditions (I) and (III) one demonstrates that the following 
inequality is met: 7 = uk/Eq ^ 1 where 7 is a so-called Keldysh parameter, and 
K = (2|£:i|)^/^. The limit 7^1 defines a weak field regime (see, e.g. [3]). Thus, in this 
paper we investigate the infiuence of the induced dipole moment D(t) on the angular 
and spectral distributions of photoelectrons in the process of single-photon ionization 
(or detachment) occurring in the weak field regime. Earlier the role of D(t) was studied 
in the process of many-photon detachment (i.e., when u <^ \ei\) from atomic negative 
ions. It was done in the weak-field [1] and in the strong-field [2] limits. 

Let us note, that the weak field condition (II) does not immediately imply the 
applicability of the perturbative approach with respect to U{r,t). Indeed, from (2)- (4) 
follows that \U{r,t)\ ~ Eo\a{Lj)\/r'^ exceeds Eq by a factor |a(ci;)|/r^. For a strongly 
polarizable target this factor can be large enough in a wide range of radial distances 
r > R, and this might lead to a non-perturbative treatment of the action of W(r, t). 
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To carry out a quantitative analysis we first postulate that under certain conditions 
the field (4) can be considered as a time-dependent perturbation which modifies the 
wavefunction of the escaping electron. The criterion of applicability of this approach 
will be formulated in the course of calculations. 

The matrix element Ai, which describes a dipole-photon transition of the electron 
from the initial bound state '(/'^(r, t) = 'i/'j(r)e~'^'* to the final state i^/ir, t) = 'i/j^\r)e~^^^ 
with the asymptotic momentum p and energy e = can be written in the following 
form: 

M = Mo + Mi + M2 + ... (7) 

The right-hand side of this equation represents the power series in U{r, t). The structure 
of the terms of the series is illustrated by diagrams in figure 1. 

nCO n n 




\u(r,t) U(r,t)\ \u(r,t) 



Figure 1. Diagrammatical representation of a single-photon transition amplitude in 
a form of perturbation series in ZY(r, t) (indicated with the dashed vertical lines). The 
double line represents the initial bound state. The solid line stands for the escaping 
electron whose wavefunction is "dressed" with all static fields. 



The zeroth order term A^o describes the direct process of photoionization standing 
for the matrix element of the operator Eorexp(— icijt)/2: 

1 r°° 

Mo = - dt 7/;* (r, t) Eor e""^* z^,(r, t) = 2716(6 - e, - u) Mq, (8) 

^ J -oo 

where 

Mo = \j drV;(-)*(r)(Eor)^A,(r). (9) 

We assume that the wavefunction iIj'^\y) (whose asymptotic form is 'the plane wave 
-|- the incoming wave') corresponds to the state "dressed" with all static fields. The 
delta- function on the right-hand side of (8) expresses the energy conservation law (6). 

The terms M.n with n>l are the corrections to the transition amplitude due to the 
n-times interaction of the photoelectron with W(r, t). Since this field explicitly depends 
on t (see (4)) then in each act of the interaction the electron energy is changed (increased 
or decreased) by u. As a result, not all terms on the right-hand side of (7) represent the 
correction to the direct amplitude M.Q., which implies the validity of (6). For example, 
evaluating M.i (the second diagram in figure 1) one finds, that it contains the terms 
proportional to 5{e — e\) and to 5{e — e^ — 2u;), which confiict with the conservation law 
(6). Therefore, the term }A\ must be ignored when constructing the correction to A^o- 

The lowest-order correction originates from the amplitude M.2 (the last diagram 
in figure 1), which contains the terms proportional to b{e — Si — lS). Applying the 
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standard technique of the time-dependent perturbation theory, one derives the following 
expression for this contribution: 

-M2|,=,,+^ = 27r5(£ -Si- u;)M2 (10) 

where 

drG,(ri,r)(Eor)V^,(r). (11) 



Here GE(r',r) is the Green's function of the photoelectron and Bq = Eq/Eq is the unit 
vector of the laser field polarization. 

Hence, within the second second order perturbation theory in U{r,t), the 
photoionization amplitude is written as: 



27r5{e -Si-u) (^Mq + (12) 

Using the amplitude (12) one can write the following expression for the spectra- 
angular distribution: 



do" p 



M0 + M2 



an 



K{u;,9)). (13) 



(14) 



Here 

stands for the differential cross section calculated without the influence of the induced 
dipole moment. Within the dipole-photon approximation the differential cross section 
dao/dQ can be presented in the well-known general form (see, e.g., [6]): 

(15, 

where ctq = o'qIuj) is the total cross section of photoionization, 6 stands for the 
angle between p and eg, P2(cos^) is the Legendre polynomial of the second order. 
The quantity (3{uj) stands for the angular asymmetry parameter [7]. The factor 
1 + l3{uj)P2{cos6) defines the profile of the angular distribution. 

Equation (13) indicates, that the profile is modified by the factor l + /i (u;, 6), when 
the field (4) is taken into account. The correction term K{u, 9) is defined as follows: 

Mo 

i^(a;,^) =2Re— . (16) 
^ ^ Mo 

Formula (13) was obtained by means of the perturbation theory. Thus, it is implied 
that the absolute value of the correction term satisfies the condition |i^'(ci;,^)| < 1. If 
otherwise, the field W(r, t) must be accounted for in all orders of the perturbation series 
(7). 

It follows from (11), that K{u!,9) is proportional to the squared induced dipole 
moment, K oc |dp = \a{uj)\'^ Eq. Therefore, for a target of a small polarizability, 
one has to apply very intensive laser fields in order to achieve a noticeable correction 
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to the photoelectron angular distribution, i.e. |i^(co',^)| ~ 1. On the contrary, for a 
strongly polarizable target (e.g., a metallic cluster or/and its anion) a strong effect can 
be expected for laser fields of a moderate strength. In section 3 we demonstrate, that 
dramatic change of the angular distribution profile can occur for Eq ~ 10^ V/cm. 

It is important to note, that the correction term (16) depends on the ratio of 
the amplitudes M2 and Mq. Thus, one can expect that it is less sensitive to the 
approximation chosen to describe the photoelectron wavefunction than the cross section 
dao/dQ. Therefore, to estimate the influence of the target polarization one can calculate 
K{uj, 6) within the framework of a comparatively simple approximation, which is 
described below. 



2.3. Analytic expression for K{uJ, 6) 

In this section we derive an approximate expression for the correction term K{uj, 6) 
using a simple, although realistic and physically clear, model. 

To start with let us assume, that the ionized target is a spherically-symmetric anion. 
Then, firstly, the ground state wavefunction ijji{r) does not depend on the angles of the 
position vector: i/ji{r) = i/ji{r). Secondly, the wavefunction 'ijj^\r) of the photoelectron 
does not contain the contribution of a long-range Coulomb field of the residue. These 
assumptions, being not of a principal nature, allow one to simplify the intermediate 
algebra. 

Making use of the partial wave expansion 

= -J2^'^~'''^^^lr.^inr)Y,Un,.)xAr), (17) 

^ Im 

where Yimin) are the spherical harmonics, 6i{p) are the scattering phaseshift, and Xpi{^) 
are the radial wavef unctions, one writes the amplitude Mq in the following form: 

Tp poo 

Mo= -2m^e'''^^U dr xAr) Ur) ■ (18) 
P Jo 1=1 

Another simplification is based on the approximation that the field U{r, t) acts on 

the escaping electron only outside the target, i.e. at the distances r > R. Additionally, 

one can assume, that the wavefunction ipi{r) is concentrated in the interior of the target, 

i.e. at r < R. Therefore, when evaluating the right-hand side of (11) one takes into 

account that ri,r2 > R > r. This allows one to relate M2 to the amplitude Mq. To do 

this let us approximate the Green's function Ge(ri,r) as follows: 

G,(ri,r) = 2ipJ2Gei{ri,r)YUn,)Y;jn) ^ ^ ^pKn)Xpz(r) 

Im Im 

Here r^/r^ is the largest /smallest of ri,r, Xpi{r) stands for the irregular solutions of 
the radial Schrodinger equation for the energy e = p'^/2. The quantity gei{fi,r) = 
Xpiij>)Xpiij<) /'P^''^>'^< is the exact radial Green's function corresponding to the 
multipolarity /. The first equality in (19) states the partial- wave expansion of the 
Green's function. The approximate equality is based on the assumption ri > r. 
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Using (19) in (11) one notices that the radial integral over r reduces to that from 
(18). Therefore, the amplitude M2 acquires the form: 

M2^-^{uj,e)Mo, (20) 

with 

-i(5i(p) 

^(00,6) = \a{uj)\'E',^— COS-' e 

dridr2 V^(-)*(r2)^ f Ge.(r2, ri) + €,^^2, ri)) ^-^^X,M). (21) 

Here the integration is carried out over the spatial region ri,r2 > R. 

The correction term (16) is related to '^{00,6) as K{uj,6) = — 2Re7(u;, ^). 

Let us comment on the evaluation of the right-hand side of (21). Introducing the 
expansion (17) and representing the Green's functions G'£j(r2,ri) and G'£_|_^(r2, ri) in 
terms of the partial-wave series similar to (19), one carries out the angular integration 
by means of the standard methods (see, e.g., [8]) and arrives at: 

K{uj,e) = -2Re7(cj,^) 



-l\a{uj)\'E',Yl 



A / A . ^ ^A A ^3(C0SI 

On Oi'o + -0i'2 - -Ol3di'2- 



5 5 cos^ 



ImBu,. (22) 



The Kronecker symbols Sij, which reduce the allowed values of the orbital momenta to 
Z = 1, 3 and /' = 0, 2, reflect the dipole selection rules applicable to each of the vertices 
in the last diagram in figure 1. The ratio P3(cos^)/ cos6' reduces to (5P2(cos6') — 2)/3 
with the help of the recurrence relation for the Legendre polynomials. 
The quantities Bw stand for the following radial integrals: 

Bw = / / dridr2 { 9eii'{r2,ri) + gs+u;i'{r2,ri) ] . (23) 

P Jr Jr n r2 \ J 

For any realistic target the exact evaluation of these integrals can be done by numerical 

means only, implying the computation of the regular Xpi{^) ^he irregular Xpi{r) 

functions as well as of the Green's functions 5'£i/'(r2, ri) and ge+oji'{f2,ri). To carry 

out the approximate integration one can substitute these quantities with the leading 

terms of their asymptotic representations. For Xpii"^) ^^^d Xpi{r) this results in 

Xpii''") ~ sin(pr — 7r//2 -|- Si{p)) and Xpi(r) ~ —e^^P^+^'>-^^\ The leading terms in the 

expansions of the radial Green's functions are: 

i &P'> SAP)&^'^< - e-'^'-< 
9e+uv{r2, ri) ~ -- , 

^ r . ^ (24) 

9eiv[r2,ri) . 

Pi r> r< 

Here P = \/2{e + uj), pi = \J —2e\, Sii{P) and 5'//(ipi) are the elements of the scattering 
matrix. For the real argument P this quantity reads as Sii{P) = 

By using the asymptotic formulae one carries out analytical evaluation of the 
integrals from (23). The final result (the leading term) for ImBui is as follows: 

lT^Bu>- (^2 _ 1)2 
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i=l,3 

;'=o,2 



where X = pR, ^ = P/p and = Pi/p. It follows from the approximations formulated 
above, that formula (25) is valid if the photoelectron momentum p and the ground state 
energy Ei satisfy the conditions: 



pR > 1, V^iR > 1. (26) 

Letting Si{p) = Si'{P) = in (25) one finds the plane- wave Born limit of ImBui. 

The presence of oscillatory terms on the right-hand side of (25) is physically clear. 
Indeed, in the limit pR > 1 the wavelength of the photoelectron is less than the size of 
the target, which is supposed to have a well-pronounced edge (e.g., the edge of an ionic 
core in a cluster). Thus, the electron experiences a diffraction at the edge when leaving 
the core. This diffraction leads to the oscillatory character of the cross section. This 
phenomenon was discussed previously in connection with various processes involving 
metallic clusters and fullerenes: the radiative electron capture by metallic clusters [9], 
the electron scattering [10], the polarizational bremsstrahlung [11]. 

Carrying out obvious summations in (22), one writes correction factor as follows: 

Kito, 9) = El [k^{uj) + K^iuj) P2(cos ^)) , (27) 

where 

1 / 4 4 \ 2 

Ko(cu) = --|a(cu)plm ( Pio + -Pi2 + -P32 1 , K2(cu) = -|a(cj)|2lmP32. (28) 

Formula (27) stresses the proportionality of K{uj, 9) to the laser field intensity 
(which is proportional to Eq) and exhibits explicitly its dependence on the emission 
angle 9, which enters, as well as in (15), via the second order Legendre polynomial. 

As a function of 9 the modulus of K{(jj^ 9) attains its maximum either at 6' = 0° 
(and 9 = 180°) or at 6' = 90°, depending on the sign of the product kq{uj)k2{u!) and 
on the relative magnitudes of k,o{^^) and K2(ci^). For further reference let us define the 
quantity 

|/«(t^)|max = ^o(t^) + «2(t^) P2(cOSI 

\ko{uj) + K2{uj)\ if Ko{uj)k2{uj) > 0, 

|Ko(ci;)| + |k2(ci;)|/2 if hq{lj)k2{uj) < and 4|Ko(ti;)| > |k2(ci;)|, (29) 
|/t2(ci;)| — |ft;o(i-^)| if f^o{^)f^2{^^) < and 4|Ko(ct;)| < |k2(i-o')|. 

which allows one to estimate the laser field strength needed to achieve a noticeable 
correction to the angular distribution (see the discussion in section 3). 

Substituting (27) into (13), taking into account (15) and expressing P2^(cos^) via 
P2(cos^) and P4(cos 9), one derives the following formula for the differential cross section 
within the lowest order of perturbation theory in U{r,t): 
da ctq 



dQ An 



(ao{uj) + a2{uj)P2{cos9) + ai{iu)P4{cos9)y (30) 
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The coefficients ao,2,4 are related to the asymmetry parameter (3 and to kq 2 via 



ao = 1 + ^oM + ^ ) , a2 = (3 + E-^ 



H2+ [fio + —j P 



The right-hand side of (30) suggests that not only the profile of the angular 
distribution can be changed, but the magnitude of the cross section, integrated over 
the emission angles, is scaled by the factor oq which can be greater or lower than one 
depending on the signs and magnitudes of /5 and /to,2- 

If the action of U{r,t) is ignored then — > 1, oq — > /3 and — ^ 0, and equation 
(30) reduces to (15). 

3. Numerical results 

Numerical analysis of the influence of the induced dipole moment on the profile of 
the angular distribution was carried out for the process of photodetachment of sodium 
cluster anions Na^. The results presented below in this section refer to = 8, 20, 40, 58 
and 92. The corresponding neutral clusters are spherically-symmetric which makes 
applicable the theory described in section 2. 

The parameters of sodium clusters and their anions, used in our numerical 
calculations, are summarized in table 1. Let us comment on these data. 

The ionic core radius was calculated using the standard relation R = VgN^''^ with 
the Wigner-Seitz radius set to 4 a.u. The values of electron affinities /a are taken 
from a recent work by Kostko [5] which contains, in particular, an excellent collection of 
reference data on various metal clusters. The other data presented in the table refer to 
the polarizability of the core and the parameters of the plasmon peak. We considered the 
dynamic dipole polarizability a{uj) of the clusters within the framework of the resonance 
plasmon approximation (see, e.g., [11, 12]): 

«(^) = «(0) 2 i • r - ^32) 

Here a(0) is the static polarizability, Ug is the surface plasmon energy, and T is the 
linewidth of the plasmon resonance. 

Due to the spill-out effect the static polarizability of a metallic cluster of a radius 
R exceeds the classical value ad = R^, which characterizes a metallic ball of the same 
radius. In the table we present the ratio a(0)/aci calculated within the RPA [13] and 
the LDA [14]) as well as the experimentally measured values [15]. The Ug values were 
obtained with the help of the relation proposed in [16] to account for the size-dependence 
of the plasmon energy: Ug = ix'Mie(l — 1.56 /R), where lvmib = 3.27 eV is the Mie plasmon 
energy in a sodium metallic sphere, and 6 = 0.54 A takes into account the spill-out 
effect. 

The size dependence of the width F of the plasmon resonance for neutral sodium 
clusters with up to several hundreds was analyzed in [17] within the RPA scheme and 
in [18] within a semiclassical approach. These two methods, being consistent for A^ > 58, 
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Table 1. Data on several spherically-symmetric sodium clusters: the ionic core radius 
R; the electron affinity /a [5], ad — R"^ is the classical static polarizability. a{0) stands 
for calculated [13, 14] and measured [15] static polarizability; LOg is the surface plasmon 
energy calculated according to [16], T is the width of the plasmon resonance peak [18]. 





R 


/a 




a 


(0)/aci 






r 




(a.u.) 


(eV) 


(a.u.) 


[13] 


[14] 


[15] 


(eV) 


(eV) 


Nag 


8.00 


0.915 


512 


1.48 


1.44 


1.72 


2.64 


0.2 


Na2o 


10.86 


1.337 


1280 


1.41 


1.37 


1.67 


2.81 


0.4 


Na4o 


13.68 


1.509 


2560 


1.10 


1.32 


1.60 


2.90 


0.7 


Nass 


15.48 


1.719 


3712 


1.24 


1.22 




2.95 


0.4 


Na92 


18.06 


1.847 


5888 


1.22 


1.20 




2.99 


0.2 



deviate noticeably for lower values. The RPA calculations reproduce the 1/R scaling 
law which, is inadequate when compared to the experimental data for < 58 (see the 
discussion in [17]). The semiclassical calculations [18] seem to be consistent with the 
experimental data in the whole region = 10 . . . 100. The values of F, presented in the 
table, were deduced from figure 1 in [18]. 

It is seen from the table, that for all A^ the plasmon energy exceeds the electron 
affinity. This fact justifies the choice of cluster anions in a view of the present study. 
Indeed, the inequality uj^ > /a leads to the additional enhancement of K{uj, 9) owing 
to the pronounced increase of the modulus of the dynamic polarizability in the region 

To conclude the discussion on the parameters of the chosen targets let us note, that 
both criteria from (26) can be easily met for Na^. The first inequality fails only in the 
narrow near-threshold region — <^ /a but is well fulfilled for a; ~ The validity 
of the second condition one proves by introducing Ei = —la.- 

The results of calculation of the quantities relevant to the correction term K{uj, 6) 
are presented in figures 2 and 3. For each anion the calculations were performed for 
the photon energies above the ionization threshold Ja and in the vicinity of the plasmon 
peak. The latter manifests itself as a resonance in the dependence of |a(u;)|^ on u, as it 
is seen from the upper panel in each graph. The dynamic polarizabilities were calculated 
via (32) with the parameters Ug and F as indicated in table 1. Given the differences 
between the calculated and experimentally measured values of the static polarizabilities 
a{0), as well as the absence of the experimental data for Nasg and Na92, for the ratio 
a{0) /li^ we used the values close to the theoretical data (see the figures captions). 

The dependences of the coefficients k,o{^) and K2(cij) on u are presented in the 
middle panels. They were computed from (28) and (25). To calculate the right-hand 
side of the latter formula beyond the plane-wave Born approximation one must know 
the scattering phaseshifts Si{p) and Si{P). To determine the phaseshifts we solved 
numerically the phase equation (see, e.g., [19, 20]) for an electron moving in a local 
potential. The model potential, which was used to describe the interaction between 
the photoelectron and the neutral residue, accounted for the following three terms: (a) 
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Figure 2. Calculation of the correction term for Na^ and for Najp as indicated. 
In each graph the upper panel represents the dependence of |Q;(a;)p on w; the middle 
panel - dependences ko,2(w) and |Av(a;)|max, eqs. (28)-(29); the lower panel - dependence 
of £'ociit('-^)! cq. (33). Parameters ujs and T arc as indicated in tabic 1. The data refer 
to the a{0)/R^ ratio equal to 1.45 for Nag and to 1.30 for Na4o . 



the short-range potential due to the Coulomb interaction with the core's electrons and 
the ionic jellium core, (b) the exchange-correlation term, which was treated within the 
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Figure 3. Same as in figure 2 but for Nagg and Nag 
a{0)/R^ is set to 1.20. 



For botli ciusters tlie ratio 



Gunnarsson-Lundqvist approximation [21], (c) the polarization potential was considered 
in the form Upoi{r) = — a(0)/2(r^+i?^)^ for r > R and Upoi{r) = if otherwise. Inclusion 
of the polarization term is important for an adequate description the low-energy electron 
scattering from metallic clusters (e.g., [22]). A more rigorous treatment accounts for 
the many-body correlations, which may be considered, for example, within the Dyson 
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equation scheme [23]. However, for the purpose of our paper it is sufficient to treat the 
polarization potential in the simple form described above. 

Comparing the middle panels in the graphs from figures 2 and 3 one notices that 
the behaviour of ^0,2(1^) (and of |K(a;)|max as well) exhibits some general trends along 
with the specific features determined by the dynamic and geometrical properties of 
an individual target. Common to all clusters is the existence of extrema for Ko,2('^) 
at uj ^ Us, which is due to the factor |a(ci;)p, see (28). The extremum can be clearly 
defined (see the curves k,o{uj) for Nag and K,2{i^) for Na92 which have distinct maxima) or 
less pronounced (both curves for Nasg) or barely seen (the case of Na4o). The variation 
in the shape of the extremum, as well as the behaviour of the curves in the regions below 
and above the plasmon resonance, can be understood by analyzing the u;- dependence 
of the right-hand sides in (28) and (25). The non-oscillatory part of the latter depends 
on the photon energy mainly through the factor (pi?)~^ = AR~'^{uj — /a)~^. This leads 
to ^0,2(1^) |a(u;)p/(u; — Ja)^. Therefore, a symmetric resonance profile of |q!(ci;)P is 
distorted due to the factor {uj — /a)~^. For a powerful and narrow resonance (as in 
the case of Na92) the distortion is comparatively small, whereas for a wide resonance 
the disproportion between the low- and the high-u; shoulders may completely change 
the behaviour of the curve in the vicinity of the resonance, as it happens for Na4o. In 
either case the factor {u — /a)~^ results in a faster decrease of |Ko,2(t^)| with u beyond 
the resonance. In the region I^^ < u < Us the curves ^0,2(1^), being enhanced by the 
factor, exhibit a non- monotonous oscillatory character (most clearly seen in Na4o, Nasg 
and Na92 graphs). This is due to the second and the third terms on the right-hand side 
of (25), which contain the oscillating factors exp(ipi?) and exp(iPi?), dependent on the 
cluster size, and the factors Si{p) = exp(i5;(p)), Si'{P) = exp {i6i{P)), dependent on the 
phaseshifts. 

The solid curves in the middle panels of each graph represent the dependences 
|'^(i-^)|max defined in (29). For a given photon energy the quantity Eq |K(cij)|max stands 
for the maximum value of \K{uj,6)\ within the interval 6 = [0°,180°]. It was already 
mentioned, that the correction term must satisfy the condition \K{llj,6)\ < 1. Thus, 
taking into account that \K{uj,6)\ < Eq\k{u)\, one can introduce a critical value 
-£^00111(1-^), which defines the upper boundary for the laser field strength consistent with 
the perturbative approach adopted in this paper: 

£'0crit(t^) = (|'t(t^)|max)~^^^ • (33) 

The corresponding intensity of the laser field (in W/cm^) can be calculated from the 
relation /ocrit ~ 3.3 x lO^^E^^j.^^, where -Eocrit is measured the atomic units. 

In the limit Eq <C ii^ocrit('^) a strong inequality \K{uj,9)\ <C 1 is valid, so that the 
photoionization process is not affected by the induced dipole moment. The infiuence of 
the latter can become pronounced for Eq < i?ocrit(i-^) resulting in \K{u!,9) \ ~ 1. 

Dependences -Eocrit(i^) are plotted in the lower panels in figures 2 and 3. For all 
clusters this function attains its minimum in the point u ^ in the vicinity of the 
plasmon resonance. The values of u, as well as of other quantities calculated in this 
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Table 2. The minima u) of the functions £'ocrit('^) (bottom panels in the graphs in 
figures 2 and 3), the values of ko,2 and l^lmax in these points, the critical strength of 
the laser field £'ocrit(w) = (|K('I')|max) and the corresponding intensity /ocrit(w), 
the ponderomotive energy Eq^^^^{uj) / Auj'^ ; is the internal anionic field. 











|t(il')|max 


^'Ocrit('^) 


^Ocrit (i^) 


E, 






(eV) 


(10«) 


(10«) 


(10^) 


(10-3 a.u.) 


(10"W/cm2) 


{IQ-^ a.u.) 


(10^3 


Nag 


2.64 


-0.470 


-0.031 


0.501 


1.41 


0.66 


17.5 


1.44 


Na2o 


2.73 


0.021 


0.022 


0.043 


4.82 


7.67 


30.8 


15.7 


Na4o 


2.72 


0.022 


-0.100 


0.078 


3.58 


4.23 


37.0 


8.64 


Nagg 


2.89 


0.151 


-0.160 


0.231 


2.08 


1.43 


44.9 


2.61 


Nag2 


2.98 


0.128 


-1.589 


1.461 


0.83 


0.23 


50.1 


0.39 



point, are presented in table 2. It is instructive to compare the critical field i?ocrit(i^) 
with the strength E'a of the internal anionic field. The latter can be related to the 
ionization potential /a through E^, = (2/a)^/^. It is seen that _Ea exceeds £'ocrit('^) by the 
order of magnitude. The strong inequality ^ _Eocrit('^), being in accordance with the 
first condition from (5), justifies the approximation adopted in this paper to neglect the 
modification of the ground state orbitals due to the laser field. The last column in the 
table contains the values of ponderomotive energy of the photoelectron in the critical 
laser field. Comparing the ponderomotive shifts and the kinetic energies e = u) — /a ~ 1 
eV, one notices that the latter is larger by several orders of magnitude. Therefore, the 
second condition in (5) is also well fulfilled. 

Hence, one can expect, that for the intensities Iq < /ocrit('^) ~ 10^^ W/cm^ the 
modification of the spectral-angular distribution due to the action of the field of the 
induced dipole moment can be quite pronounced. To carry out the corresponding 
quantitative estimates one can use the formulae obtained within the framework of 
perturbation theory. 

The calculated profiles of the angular distributions are presented in figure 4. In 
the present paper we have not made an attempt to compute the angular asymmetry 
parameter i3{uj) for sodium anion clusters (neither have we found its values in the 
literature). Instead, for each target we carried out the calculations for the three values, 
/3 = —1, 0, 2, aiming to demonstrate that the modification of the profile takes place for 
any (3 from the allowable interval [—1, 2]. 

In each graph from figure 4 the thick solid curve represents the unperturbed profile 

''^■sin^^ for/3 = -l. 



F^(e) = l+/3P2(cos( 



2 

1 ioi (5 = 0, (34) 

3 cos^ e for /3 = 2. 

The emission angle 9 is measured with respect to the vector Eq, which is indicated in 
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Figure 4. Profiles J-'i3{0) of tiie angular distribution (35) calculated for different 
sodium cluster anions and for the three values of asymmetry parameter (3 (as 
indicated), and for several values of the ratio r] = Iq/Iq ait (see the legend in the 
graph /3 = 2 for Na^). The data refer to the photon energies ui indicated in table 2. 
Thick black curves in each graph represent the field free distribution Fp{6), eq. (34). 
The length of line segment connecting the origin and a given point on the curves is 
equal to the value oi J-'plO) or Fp{9) (in absolute units) in the corresponding direction. 
The emission angle 9 is measured with respect to Eq, which is shown in the graphs 
with (3 = —1. In each graph the spacing between ticks is equal to 1 for both directions. 

the graphs with (3 = —1. Other curves correspond to the modified profiles 

rp{e)^Fp{e)(i + ElU + K,p,{cose))) =F,ie) A +^ ^o + /^2P2(cosg) \ ^^^^ 

^ ^ ^ ' \ 1^1 max / 
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(see (13), (27) and (33)) calculated for several values of the ratio rj = /o//ocrit < 1- For 
each anion the data refer to the photon energy u. 

Figure 4 demonstrates, that even for comparatively low intensities (corresponding 
to ?7 > 10~^) the changes in the profiles are already visible. For t] < 1 they become much 
more pronounced. It is seen that variation of the profile depends, both qualitatively and 
quantitatively, on the value of f3 and on the choice of the target anion which enters via 
the coefficients kq and k,2- 

In the case jS = —1 the field-free profile Fp{6) is proportional to sin^ 6' yielding 
-F_i(0°) = as the minimum and F_i(90°) = 1.5 as the maximum values within the 
interval 6 = [0°, 90°]. The modified profile can be written as jF_i(6') oc sin^ 6 (a + 6cos^ 6) 
with a = 1 + 77(2^0 — '^2)/2|K|max and h = 377^2/2 |/t|jnax- One easily shows, that for 
\a\ < \b\ the function T-i{9) acquires additional extremum. However, this inequality 
cannot be met for the considered targets (see the values of kq, 1^2 and l^lmax in table 
2) and for r/ < 1. Therefore, the pattern of the modified angular distribution is similar 
to that obtained without the infiuence of the induced dipole. Namely, the strongest 
emission occurs in the directions perpendicular to the field Eq and there is no emission 
for 6' = 0° and 180°. On the whole, the distortion increases with 77, attaining the largest 
values at ^ = 90°. The trend of the largest distortion depends on the sign of n^. if 
> (as in the cases of all anions but Nag) then jF_i > -F-i whereas negative kq leads 
to a shrinking of the modified profiles (as it happens for Nag^). 

For /5 = the unperturbed angular distribution is isotropic. The modified profile 
reads as jFo(^) oc a + 6 cos^ 9, with a and h defined as above. The way the profile evolves 
with the increase of 77 to a great extent depends on the sign of H2. For ^2 > the emission 
along ±Eo dominates over the emission in the perpendicular directions (see the case of 
Na2o), whereas for negative K2 the perpendicular direction becomes preferable. The 
graphs from the central column in figure 4 demonstrate, that for 77 ^ 1 the distortion 
can become dramatic both in form and in magnitude. 

In the limit (3 = 2 the modified profile behaves as T2{0) oc cos^ 9 {a + bcos^ 9). The 
additional extrema appear if the following two conditions, ab < and \a\ < 2\b\, are 
met. These conditions can be satisfied in the limit rj ^ 1 for all considered anions but 
Na2o- Most explicitly it is seen in Najg and Na^2 graphs, where the curves with rj = 1 
represent the four-petal angular distributions with the maximum emission at ^ = 45° 
and 135°. The same pattern characterizes Nag , although in this case it is hardly visible 
due to a strong decrease in the magnitude of ^2(6*) for 77 = 1. 

Finally, let us briefiy discuss the influence of the induced dipole moment on the 
cross section integrated over the emission angles. Recalling (30) and (31), one finds 
the following formula for the ratio of the cross sections with, o'{uj), and without, o'q{uj), 
account for the correction: 



Figure 5 presents the ratios o-{uj) / cro{u) versus to. For each anion the ratio was calculated 
using the dependences /to(co') and ko(u;), presented in figures 2 and 3, and for three 




(36) 
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Figure 5. Ratio (7{uj)/ao{uj) (sec (36)) calculated for different sodium cluster anions 
and for the three values of asymmetry parameter /3 (as indicated in the top left graph). 
The data refer to the laser field intensity Iq = 2.1 x 10^" W/cm^. 



values of the asymmetry parameter (as indicated). The intensity of the laser field was 
set to Jo = 2.1 X 10^° W/cm^, which corresponds to the field strength £"0 = 8 x 10~^ 
a.u. Comparing the latter value to the critical field Eocriti^) (see figures 2 and 3) one 
notices, that for all targets and in the whole range of u the parameter 77 = E^/ Eq^j.^^{uj) 
is less than one. Thus, formula (36) is applicable for the quantitative estimate of the 
correction due to the induced moment. 

As one can see from figure 5, the largest deviation (up to ~ 30%) of cr{uj) from 
ao{uj) occurs for Na^ and Nagg in the vicinities of the plasmon resonances, where r) 
reaches its maximum values: rj ^ 0.57 for Nag and rj ^ 0.95 for Na^2- ^^e case of 
Na92 a strong inequality |k2| ^ \i^o\ is valid (see figure 3 or/and table 2). Therefore, 
the magnitude of the correction term from the right-hand side of (36) becomes sensitive 
to the choice of the asymmetry parameter. As a result, the ratio a{uj)/ao{uj) strongly 
varies with (3. On the contrary, for Na^ the relation \k2\ -C \ko\ is valid, leading to a 
weak dependence of the ratio on p. 

For the anions Najg and Nagg the largest values of rj equal to 0.05 and 0.15, 
correspondingly. Hence, the intensity Jq is small compared to /ocrit('^) oc Eq^j.^^{lj). 
As a result, the deviation of cr{uj) from the field-free value cro{uj) is less pronounced. 
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4. Summary 

In summary, we have demonstrated that the dipole moment, induced at the residue 
by the incoming laser field, can strongly infiuence the photoionization process. The 
additional time-dependent long-range field due to the induced moment, acting on 
the escaping electron, modifies the spectral-angular and spectral distributions of 
photoelectrons. 

The effect itself, as well as its quantitative treatment, is quite sensitive to the 
choice of the target and to the parameters of the laser field. We have demonstrated, 
that for a strongly polarizable target (e.g., metallic cluster anions) strong modifications 
of the characteristics of a single-photon ionization process in the vicinity of the plasmon 
resonance (typically, Us ~ 2 — 3 eV) can be achieved by applying laser fields of 
comparatively low intensities /q ~ 10^° — 10^^ W/cm^. 

From a theoretical viewpoint, the advantage of this regime is that the field of 
induced dipole moment can be treated perturbatively, whereas the infiuence of the laser 
field on the photoelectron (a ponderomotive effect) and on the ground state of the 
target can be ignored. We have shown, that even in this weak-field regime one can 
expect dramatic changes in the profile of angular distribution of photoelectrons as well 
as in the u dependence of the photoionization cross section. To a great extent, these 
changes depend on the parameters of the target, like its radius, the static and dynamic 
polarizabilities, the energy of the plasmon resonance peak and the affinity (or, the 
ionization potential in the case of a neutral target). Additional variety can be obtained 
by changing the intensity and the frequency of the laser field. 

The analysis, which we have carried out, is based on simple models adopted for the 
description of the photodetachment process (the plasmon resonance approximation) and 
of the interaction of the escaping electron with the static field of the residue (via the LDA 
-|- long-range polarization potential). We have also restricted ourselves to the case of 
spherically-symmetric targets. A more accurate quantitative treatment must include 
more sophisticated approaches, which take into account the many-body correlation 
effects intrinsic to both the photodetachment and the electron escape. However, these 
will leave untouched the physics behind the basic phenomenon discussed in our paper. 

It is possible to study the predicted phenomena by means of modern experimental 
techniques. The latter include the production and manipulation of the beams of 
cluster anions of various types and sizes, the energy and angle resolved photoelectron 
spectroscopy, the availability of the laser fields of various intensity, frequency and pulse 
duration [4, 5] . We believe, that with all these components at hand, soon it will become 
possible to test the theory against experimental data. 
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